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ABSTRACT 

Atomic diffusion may lead to heavy element accumulation inside stars in certain specific 
layers. Iron accumulation in the Z-bump opacity region has been invoked by several authors to 
quantitatively account for abundance anomalies observed in some stars, or to account for stellar 
oscillations through the induced K-mechanism. These authors however never took into account 
the fact that such an accumulation creates an inverse /i-gradient, unstable for thermohaline 
convection. Here, we present results for A-F stars, where abundance variations are computed 
with and without this process. We show that iron accumulation is still present when thermohaline 
convection is taken into account, but much reduced compared to when this physical process is 
neglected. The consequences of thermohaline convection for A-type stars as well as for other 
types of stars are presented. 

Subject headings: stars:abundances, stars:oscillations, hydrodynamics, convection, asteroseismology 



Introduction 



Th e importance of atomic diffusion (jMichaud 
I1970I) inside stars is now well established: not 
only can it modify the atmospheric abundances, 
as observed in the so-called "chemically peculiar 
stars" , but it can also have strong implications 
for the stellar internal structure. The variations 
with depth of the radiative accelerations on spe- 
cific elements can lead to their accumulation or 
depletion in various layers inside the stars. In 
this framework, the accumulation of iron (together 
with less important metals like nickel) may have 
strong consequences, which are different according 



to the stellar type. 

In main sequence A stars, it may lead to the 
creation of an extra convective zone due to an iron- 
induced increase o f opacity (jRicher et al. l200d: 



Richard et al ] l200lh . If overshooting leads to com- 



plete mixing between this extra convective zone 
and the one due to helium, helium settling may 
be slowed down in the outer stellar layers. Such 
an effect could expl ain the existence of e volved os- 
cillating Am stars ( Turcotte et al. 2000l ). The oc- 
currence of an iron convective zone may also help 
explaining the excitation of g-modes in 7 Doradus 
stars, through the "convective blocking" process, 
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as suggested by Guzik et al. ( 2000l ). Furthermore, 
in some cases iron accumulation can trigger stellar 
pulsations through the iron-induced K-mechanism. 
This process has been invoked to explain main se- 
quenc e B-type pulsators like SPB and Cephei 



stars (IPamvatnvkh et al 



2004 iBourge fc Alecian 



icd for the case of accretion of metal-rich matter 
onto main sequence stars. This situation may 
occur for exoplanet host stars, if they accrete 
hydrogen poor matter in the process of plane- 
tary formation. Because of thermohaline convec- 
tion, and contrarily to preliminary conclusions 



20061 ). Charpinet et al. (1996) predicted that the found in early papers ( Santos et al. 2001 ). the 



extreme horizontal branch stars referred to as sdB 
stars could also be destabilized due to iron-induced 
K-mechanism. The oscillations of these stars were 



independently observed (e.g. lO'Donoghue et al 
19971 ) . The qualitative success of the iron accumu- 



lation theory is impressive. However when quanti- 
tatively comparing models and observations, some 
difficulties remain. We will further discuss these 
difficulties in section [5] 

In previous stellar modelling done in this frame- 
work, the question of the stability of a heavy iron 
layer lying above layers with smaller mean molec- 
ular weights was not addressed. This represents 
a serious problem, because of the existence of lo- 
cal inverse /i- gradients, which are highly unstable 
agai nst th ermohaline (double diffusive) convection 
(see lVauclair.2004 ') . The induced mixing has to be 



taken into account in the computations, otherwise 
the results which are obtained do not correctly 
represent the real situation found inside the stars. 

The importance of thermohaline convection has 
been invoked in a few astrophysical cases. The 
first case which was pointed out c oncerned the 
•^He burning regions in gi ant stars (Ulrich 1972t 
Charbonnel fc ZahnI 120071 ) . Situations with ^He 
enhancement in outer stellar layers were also dis- 
cussed in several studies. This enhancement could 



be due to mass transfer (|Stothers fc Simpnl 119691) 
or to helium diffusion in a stellar wind ( Vauclair 



19751 ) . The observations of main sequence helium- 
rich stars, with effective temperatures around 
20000K, show that helium is enhanced in average 
by a factor of two. This demonstrates that ther- 
mohaline mixing is important, otherwise the he- 
lium overabundance would be much larger. How- 
ever, thermohaline convection mixing cannot be 
extremely efficient since this would lead to no he- 
lium overabundance. In this case, the reason for 
this lack of efficiency may be due to the pres- 
ence of a magnetic field. The same kind of pro- 
ce sses may occur inside r oAp stars, as discussed 
bv lBalmforth et all (|200ll ). 

Thermohaline convection has also been stud- 



accreted matter does not remain inside the con- 
vecti ve zone, but diffuses downwards (jVauclair 
20041 ). Another case concerns carbon enhanced 
main-sequence stars (CEMPs), which are as- 
sumed to have suffered some accretion of ma- 
terial coming from an AGB companion. This 
heavy matter falls down inside the star due t o 



thermohaline convection ( Stancliffe et al.l 120071 ). 



However, owing to the stabilizing /^-gradient in- 
duced by helium diffusion, some of it can re- 
main in the outer layers and may account for 
the observed abundances (iThompson et al. 1 120081: 



Stancliffe fc Glebbeek|[2008t IStancliffell2009l) . 

In this paper, we ffist give a general discussion 
of thermohaline convection in stars (section[2]) and 
show how this physical process can be introduced 
in the stellar evolution computations (section [3]). 
Then results are presented for heavy elements ac- 
cumulation inside A-type stars. We mainly con- 
centrate on the evolution of a 1.7 Mq star, but 
also give results for 1.5 Mq and 1.9 Mq stars (sec- 
tion jH). For these compvitatio ns, the TGEC code 
(cf. iRichard et al. 2004; Hui- Bon-Hoal l2008f ) was 
modified to include the computations of radiative 
accelerations using the SVP (standing for Single 
Valued Parameters ) appr oximation proposed by 
Alecian fc LeBland (|2002l) . This modified version 



of the TGEC code was first used to compute stellar 
evolution model s (without introduci ng thermoha- 
line convection i Theado et al. 2009f )) which were 
compared to the results obtained by the Montreal 
group. These computations were repeated while 
including thermohaline convection below the iron 
accumulation layer. Due to the stabilizing effect 
of helium settling, which fortunately occurs in the 
same stellar region, iron can still accumulate, but 
the abundance profiles are strongly modified. We 
also discuss the potential important consequences 
of this process for other types of stars in section 
[51 A short conclusion will then follow. 
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2. Thermohaline convection 

Thcrmolialinc convection is a well-known pro- 
cess in oceanography: salted water layers lying 
above fresh water ones are rapidly mixed down- 
wards even in the presence of stabilizing tempera- 
ture gradients, due to the different diffusivities of 
heat and salt. When a warm salted blob falls down 
in cool fresh water, the heat diffuses out more 
quickly than the salt. The blob goes on falling 
due to its weight until it mixes with the surround- 
ings. This leads to the formation of so-called "salt 
fingers" . 

A similar kind of convection occurs in stel- 
lar radiative zones in the presence of positive 
/x-gradicnts, that is /^-values increasing upwards. 
An "iron layer" lying above hydrogen rich medium 
will lead to "iron fingers" , because when a blob 
falls down heat diffuses between the blob and its 
surroundings more quickly than iron ions. With 
reference to oceanography, this process is called 
"thermohaline" convection. It is also sometimes 
referred to as double diffusive convection, because 
it is directly related to the presence of two scalars 
which diffuse on very different time scales. 

Another kind of double diffusive convection oc- 
curs in the reversed case, in the presence of desta- 
bilizing temperature gradients and stabilizing /i- 
gradients like, for example, at the edge of con- 
vective cores. This process, quite different from 
thermohaline convection, is generally referred to 
as semi convection. It will not be considered here. 

In the case of positive /x-gradients, the medium 
can become dynamically unstable if (Ledoux cri- 
terion): 



^ r.rit ^ I. 



"^ad - V < 



(1) 



where cj) = {d \n p/d In fi) and S = {d In p/d 
In T). When Vcrit vanishes, marginal stability is 
achieved. For positive or null values of Vcru , the 
medium is stable against dynamical convection. 
Then thermohaline convection occurs on a time 
scale which is typically of a few thousand years, 
long compared to the dynamical time scales but 
short compared to stellar lifetimes. 

Fingers may form if the following condition is 
verified: 



with r — D^/Dt — tt/t^ where Dt and I?^ are 
the thermal and molecular diffusion cocfhcicnts 
while tt and are the corresponding time scales. 

In stars the value of the t ratio is typically on 
the order of 10~^" if one assumes that is the 
"atomic" diffusion coefficient. However, it can be- 
come larger when the shear flow instabilities which 
induce mixing between the edges of the fingers 
and the surro undings are taken into account (see 
Vauclaii|[200l . 

The effects of thermohaline convection as a mix- 
ing process in stars are far from trivial. Many 
detailed studies related to the occurence of this 
physical process in water have been p ublished. For 
example, iGargett fc Ruddickl (|2003l ) gave precise 
comparisons between numerical simulations and 
laboratory experiments. However the stellar case 
is different since mixing occurs in a compressible 
fluid. 

Two differ ent parametr isatio n recipes, given re- 
specti vely bv lUhichI (|l972[) and lKippenhahn et a 



(2) 



(|i98n[ ) can differ by up to two orders of magnitude: 
this illustrates that treating thermohaline convec- 
tion as a simple diffusion process may lead to in- 
accurate results. The basic problem here concerns 
the vertical shear flow instability which occurs be- 
tween the fingers and the inter-fingers medium. 
This instability leads to local turbulence which in- 
creases the mixing at the edge of the fingers. Con- 
sequently, a process of self destruction appears for 
the blobs, so that fingers eventually reach a regime 
where they cannot form anymore: this effect is 
taken into account in the Kippenhahn et al. pro- 
cedure. 

In the present paper, we take into account the 
fact that the thermohaline mixing time scales are 
short compared to the main sequence evolution 
time scales, and that mixing stops when the posi- 
tive /i-gradients vanish. In our computations (see 
discussion below), we model rapid mixing so that 
the stellar medium always adjusts to keep either 
fiat /i-profiles or /i- values increasing downwards. 
As helium settles downwards, it creates a stabi- 
lizing /i-gradicnt which in turn may allow some 
heavy elements accumulation. 

3. Computations 

We computed stellar models to describe the 
evolution of A-type stars while using the Toulouse- 
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Geneva Evolution Co de (TGEC). The co de is 

chard et all \2004) and 
except for the computa- 



described in d etail 
Hui-Bon-Hoal (|2008h . 



tions of opacities and atomic diffusi on which have 
unde rgone major improvements (jXheado et al 
20091) . 



3.1. Opacities and radiative accelerations 

fn the present computations, important im- 
provements were introduced in the TGEC code 
concerning the opacities and radiative accelera- 
tions. The opacities were co mputed using the 
OPCD v3.3 codes and datsQ (jSeatonl [200I . ft 
allows to compute self consistent Rosseland opac- 
ities taking into account the detailed composition 
of the chemical mixture. The opacities were recal- 
culated at each time step. 

To more realistically treat atomic diffusion 
in the TGEC code, the radiative accelerations 
on C, N, O, Ca and Fe were included foUowini 



the improved version (jLeBlanc fc AlecianI 1200 
of the semi-analyt i cal pr escription proposed by 
Alecian fc LeBland (|2002l) . This method allows 
for very fast computations of radiative accelera- 
tions with a reasonable accuracy. The SVP ap- 
proximation may be implemented in existing codes 
in a simple way. There are much less data to pro- 
cess than for detailed radiative acceleration calcu- 
lation, because complete and detailed monochro- 
matic opacites for each ion are not needed. A 
new grid of SVP-parameters, well fitted to the 
stellar mass range considered in this work, was 
computed following the procedure described in 
LeBlanc fc AlecianI (|2004 . 

The diffusion of H, He and the five metals men- 
tioned above were included in the stellar mod- 
els. The abundance of the other elements were 
assumed to be constant throughout the star. The 
radiative accelerations of these five heavy elements 
as well as the detailed abundances for all seven el- 
ements were computed at each diffusion time step, 
chosen to be about 50 times smaller than the evo- 
lution time step. 

3.2. Convection and Mixing 

For our modelling, the convective regions were 
instantaneously homogenised. The HI and HcH 



^The OPCD_3.3 packages is available on the following web- 
site: |http://cdsweb.u-strasbg.fr/topbase/op.html| 



convective regions were supposed close enough to 
be connected by overshooting and mixed together. 
On the other hand, the iron convective region 
which may appear in some models in much deeper 
regions, was supposed disconnected from the sur- 
face convective zone. 

The models were evolved from pre-main se- 
quence up to hydrogen core exhaustion. They 
were assumed homogeneous on the pre-main se- 
quence and atomic diffusion was introduced at the 
beginning of the main sequence. 

We first computed models without thermoha- 
line convection. To avoid the appearance of steep 
and unrealistic abundance gradients at the transi- 
tion between radiative and convective regions, we 
introduced mild mixing at the bottom of each con- 
vective zone. This mixing was modeled as a diffu- 
sion process through a coefficient Dturb introduced 
in the chemical transport equation: 



dcj 
dt 



j,2 Qj. 



r pDturb 



dt 



(3) 



Ci represents the mean concentration of a species 
i, p is the density, r is the radius, Dturb is given 
by: 

, r - Tczb ' 



Dturb = Dczb cxp In 2 



A 



(4) 



Dczb and rczb are respectively the value of Dtmb 
and the value of the radius at the base of the con- 
vective zone, A is the half width of the mixing 
region. Dczb and A are free parameters, choosen 
to produce a mild mixing on a small extent. The 
value of Dczb is taken equal to fO^cm^.s^^ and A 
is fixed to 0.5% of the stellar radius below the sur- 
face convective region and to 0.05% of the stellar 
radius below the iron convective zone. 

3.3. Modelling thermohaline convection 

As discussed in the previous sections, thermo- 
haline convection arises in regions where the mean 
molecular weight increases towards the surface. It 
leads to mixing on short time scales, compared 
to those of evolution, and it disappears when the 
mean molecular weight gradient vanishes. To sim- 
ulate this phenomenon, we found that the most 
efficient and simple way was to model the mixing 
as usual, provided that the chosen diffusion coef- 
ficient is large enough to flatten the /i-profile any 
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time a destabilizing /i-gradicnt is created (at each 
time step). 

Th e diffusion coefficient proposed bv lKippenhahn et al 
()l980l) is convenient in this framework: 



Dth = 



H„ l6acT^,d\n 



Vad - V CpKp^ 



dr 



(5) 



where Hp is the pressure scale height, a the radia- 
tive pressure, c the speed of light, T the tempera- 
ture. Vad the adiabatic gradient, V the real tem- 
perature gradient, cp the specific heat at constant 
pressure, k the opacity and /i the mean molecular 
weight of the material. 

We used this prescription in our computations 
and checked that the /^-gradient remained always 
null or negative at each time step in all our models. 

4. Results 

We computed models for 3 different masses: 
1.5, 1.7 and 1.9Mq. Each set of models includes a 
first series of computations without thermohaline 
convection and a second series including it. The 
initial metal mixture used in the se computations is 
the so lar mixture presented by iGrevesse fc Noeld 
(jl993l ) with the following values as H and He ini- 
tial mass fractions: Xq = 0.7112 and Yq = 0.2714. 
The initial iron mass fraction was consequently 
Xo(Fe) = 0.1149 x IQ-^. 

Figure[T]presents the evolutionary tracks for the 
3 sets of models. Inclusion of thermohaline con- 
vection does not modify these tracks in a sensible 
way: the two curves coincide. The crosses show 
the positions of the models which are discussed 
below. 

4.1. Results obtained for l.TM© models 
with and without thermohaline mix- 
ing 

For the discussion here, we specifically focus on 
the case of 1.7Mq models. Figures [2] to [6] compare 
the results obtained with (solid lines) and with- 
out (dashed lines) thermohaline mixing, at vari- 
ous epochs during the H burning phase (crosses in 
Fig. [T|) . Figure [2] displays the ratio of the iron 
mass fraction at various times to that of its initial 
value. Figure [3] shows the opacity profiles at these 
same evolution times. Figure [3] presents the differ- 
ences between the radiative and adiabatic gradi- 
ents, which allows to locate the convective regions. 



1.4 - 



1.2 - 
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3.9 3.85 
l°g T.,, 

Fig. 1. — Evolutionary tracks for three different 
masses, including atomic diffusion and thermoha- 
line convection. The crosses locate the position of 
the models shown in Figures[2]to[6l Those situated 
on the 1.7 Mq track locate models at 62 Myrs, 86 
Myrs, 299 Myrs, 403 Myrs, 588 Myrs, 788 Myrs, 
988 Myrs, 1188 Myrs and 1388 Myrs. On the 1.5 
and 1.9 Mq track, only two models are presented, 
respectively at 838 Myrs and 632 Myrs. 
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Figure [5] displays the mean molecular weight pro- 
files and Figure [6] displays the mean molecular 
weight gradient, for the interval 4.5 < logT < 6.5. 

The models reach the main sequence with 
chemically homogeneous envelopes. Then, due 
to atomic diffusion, helium sinks and creates a 
stable (negative) //-gradient below the convec- 
tive envelope, whereas the radiative acceleration 
on iron induces an increase of its abundance in 
the upper stellar layers, and more specifically in 
the iron-peak elements opacity bump centered at 
logT ~ 5.2. For the 1.7M0 case, the models with 
and without thermohaline mixing begin to diverge 
at about 400 Myrs. 

At the age of 62 Myrs, in both cases the 
iron accumulation reaches a factor of 7 below 
the surface, which does not alter the opacity, 
and a factor 4 in the opacity bump region, 
which leads locally to a small opacity increase. 
Consequently, the radiative gradient (V^ad = 

increases, but it re- 
mains smaller than the adiabatic gradient so that 
the corresponding region is still convcctively sta- 
ble. At this time (62 Myrs) the model only has 
a single surface convective zone, due to the ion- 
ization of hydrogen. As can be seen in Figure 
m from the close-to-zero value of {Vrad — ^ ad) 
around logT ~ 4.6, the Hell convective zone has 
just disappeared, due to He settling. At this point, 
the iron accumulation is not large enough to sig- 
nificantly alter the /i-profile. Here, the effect of the 
He-induced /i-gradient is indeed much more im- 
portant than the iron contribution. During stellar 
evolution, as the atomic diffusion proceeds, the re- 
gion where the helium-induced stable ^-gradient 
is the steepest slowly deepens inside the star. As 
soon as it occurs below the iron opacity bump, the 
effects of the iron accumulation on the /.(-profiles 
become important. 

At 86 Myrs, iron accumulation exceeds a fac- 
tor of 19 and affects the //-profile. It is respon- 
sible for the change in the slope of the /x-profile 
at logT ~ 5.3. The iron enrichment leads to a 
significant opacity increase and therefore to an in- 
crease of the radiative gradient which now exceeds 
the adiabatic one. In both cases, with or without 
thermohaline mixing, a new very narrow convec- 
tive region appears around logT ~ 5.2. 

At 299 Myrs, this iron convective region be- 
comes significant and leads to chemical homogc- 



nization between log T c± 5.2 and 5.4 as seen in the 
/X- and (Fe/Feo)-profiles shown. The homogeniza- 
tion of the iron convective zone leads to an expan- 
sion of the iron accumulation region and reduces 
the iron maximum value from a factor of 19 to a 
factor of 12.5. From the beginning of the main- 
sequence until 299 Myrs, the //-gradient remains 
null or negative in most models, except below the 
surface convective zone where it sometimes takes 
a slightly positive value, due to iron enrichment in 
the outer regions. 

After 400Myrs, the /(-gradient induced by He- 
settling is then located below the iron accumula- 
tion region and the effects of iron accumulation 
become visible. From now on, models including 
thermohaline convection below the iron accumula- 
tion layers significantly differ from those without 
thermohaline convection. 

During the subsequent stellar evolution, radia- 
tive levitation leads to iron accumulation in the 
iron opacity bump region. In models without ther- 
mohaline convection, this iron accumulation dras- 
tically increases with time, reaching a factor of 
20 at 400 Myrs and up to a factor of 95 at 1388 
Myrs (Fig. [21 dashed fines). The consequences on 
the opacity profiles can be seen in Figure [3l and 
those on the radiative gradients in Figure [H It 
also leads to a new convective zone which persists 
during most of the main sequence lifetime. 

Iron accumulation in these specific layers leads 
to a spectacular increase of the /i-values (Fig. [SJ 
dashed lines). After 400 Myrs, "/t-bumps" develop 
with a rapid increase slightly above logT ~ 5.2 
and a steep decrease below logT ~ 5.4, which is 
highly unstable. The corresponding /i-gradients 
are shown in Figure [HI 

In the models including thermohaline convec- 
tion induced by these unstable /t-gradients, iron 
accumulation in the opacity bump region is dras- 
tically reduced and never exceeds a factor of 15. 
It is never completely suppressed however, due to 
the stabilizing effect of helium settling. The con- 
sequences on the opacity profiles, radiative gradi- 
ents, /(-values and their gradients can be seen in 
Figure [3] to [H solid lines. The evolution of the 
internal stellar structure in the two cases studied 
here is sensibly different. 

Contrarily to the models without thermohaline 
mixing, those which include it do not have per- 
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sistent iron convective zones, although marginal 
convection appears and disappears several times 
during stellar evolution. Actually, in models in- 
cluding thermohaline mixing, the iron abundance 
in the opacity bump always stays close to the crit- 
ical value for the onset of convection. When iron 
accumulation exceeds this critical value, convec- 
tive mixing occurs and reduces the iron accumu- 
lation, making this region stable again. Note that 
when this occasional convective zone appears, it 
has a smaller extent than the iron bump itself, so 
that only part of the iron enriched zone is homog- 
enized. This is the reason why a small iron peak 
remains around logT ~ 5.2 in the 58 Myrs, 788 
Myrs and 988 Myrs models. 

4.2. Comparisons between the 1.5, 1.7 and 
I.9M0 cases 

The results obtained for 1.5 and 1.9Mq mod- 
els are qualitatively similar to those of the 1.7M0 
models. To illustrate the efficiency of thermoha- 
line convection in these models, we present in Fig- 
ure[7]the molecular weight and the iron profiles ob- 
tained with and without thermohaline mixing for 
one model for each mass, located at the middle of 
the main-sequence (cf. Fig. [T]). The radiative ac- 
celeration on iron increases with the stellar mass, 
so that, when no thermohaline convection is taken 
into account, the iron accumulation is larger in 
more massive models. The consequent "/i-bump" 
is also more prominent, as can be seen in Figure 
[7] (dashed lines). In the results obtained including 
thermohaline mixing (solid lines), we can see that 
iron accumulation is again drastically reduced. In 
the 1.9Mq model, it does not exceed a factor of 
10 whereas it goes up to 150 when thermohaline 
mixing is neglected. 

Figure [8] presents the position of the boundaries 
of the convective zones inside the models with and 
without thermohaline mixing for the three masses 
studied here. The crosses locate the boundaries of 
the iron convective zone during phases with mix- 
ing episodes. When the convective zone settles 
durably, its boundaries are represented with solid 
lines. 

For the three stellar masses considered, an iron 
convective zone resulting from a diffusion-induced 
iron accumulation rapidly appears after the be- 
ginning of the main sequence. In models without 
thermohaline mixing (Fig. [51 upper panels), a first 



phase occurs during which convective episodes are 
followed by convectively stable periods. When 
convection appears in the opacity bump region, 
the iron peak decreases due to convective mixing. 
Consequently the opacity also decreases and the 
radiative gradient becomes smaller than the adi- 
abatic gradient. The region becomes radiatively 
stable again. In the 1.5M0 models, this succes- 
sion of convective and radiative episodes goes on 
until the end of the main sequence. In the 1.7 and 
I.9M0 models, after a few hundred megayears, 
iron accumulation has become large enough (more 
than a factor of 15), so that convective mixing 
cannot reduce the iron abundance enough to stop 
convection. The iron convective zone then persists 
during the rest of the main sequence lifetime. 

In models including thermohaline convection, 
the situation is different. Here again, there is a 
first phase during which mixing episodes alternate 
with convectively stable periods. Afterwards the 
onset of thermohaline mixing strongly modifies the 
convective (un)stability of the iron opacity bump 
region. Iron accumulation is reduced and leads, 
most of the time, to a radiative gradient slightly 
below the adiabatic gradient. When, during stel- 
lar evolution, the accumulation of iron exceeds, 
somewhere in the opacity bump region, the criti- 
cal value for the onset of convection, the induced 
mixing rapidly reduces the iron abundance and the 
radiative stability is restored. As a result the iron 
convective zone does not persist in these models. 

5. Discussion 

Thermohaline convection is expected to take 
place inside stars any time inverted /i-gradients 
are built. This happens in all types of stars 
where heavy element accumulation occurs in spe- 
cific stellar layers, either due to radiative levita- 
tion or to nuclear reactions during stellar evolu- 
tion. It also occurs in case of accretion of hy- 
drogen poor material onto stellar surfaces, due to 
planet swallowing in early phases, or to accretion 
of matter from a companion. Here we presented 
how introducing thermohaline convection modi- 
fies the results obtained for the accumulation of 
iron due to diffusion inside A-type stars. Some 
other cases where this process has to be intro- 
duce d were already discussed in previous papers 
(e.g. ICharbonnel fc ZahnI (|2007t ): IStancliffe et al ' 
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Fig. 2. — Iron accumulation in two 1.7 Mq models with atomic diffusion at various evolutionary steps. 
Fe/Feo represents the ratio between the iron mass fraction and its initial value. The dashed lines are for 
models without thermohaline convection and the solid lines for models including it. 
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Fig. 3. — Opacity profiles in the two 1.7 Mq models presented in Figure [2l Dashed lines: without thermo- 
haline convection; solid lines: with thermohaline convection. 
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Fig. 4. — Difference of the radiative and adiabatic gradients in the 1.7 M© models presented in Figure [2] 
Convective regions appear when this difference is positive. Dashed hnes: without thermohahne convection; 
sohd hnes: with thermohahne convection. 
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Fig. 7. — Iron profiles (left panels) and /i-profiles (right panels) in two 1.5 and 1.9Mq main sequence 
models (with a central hydrogen mass fraction Xc ~ 0.5) including atomic diffusion. Dashed lines: without 
thermohaline convection; solid lines: with thermohaline convection. 
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Fig. 8. — Convcctivc zones in models with (lower panels) and without (upper panels) thermohalinc mixing. 
The vertical axis, log(AM/M), represents the outer mass fraction. For each plot, the upper solid lines locates 
the bottom of the surface (H-He) convective zone. During stellar evolution, a new convective zone induced 
by iron accumulation appears deeper in the stellar interior (typically around —7 < log(AM/M) < —5). In 
some models this iron convection zone settles durably while other models undergo mixing episodes followed 
by convectively stable phases. During the mixing episodes, the boundaries of the iron convective zone are 
represented by crosses. During the convectively unstable phases, they are represented by solid lines. 
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Before giving our general conclusions in 
the next section, we first propose a discussion of 
three other cases where thcrmohalinc convection 
should be introduced and for which we expect im- 
portant consequences. 

5.1. The SPB / /3 Cephei stars 

The SPB (slowly pulsating B stars) and the /? 
Cephei stars are main sequence B type pulsators. 
The SPB stars oscillate with periods between 0.5 
to 5 days, while the (3 Cephei stars possess 2 to 
8 hours-oscillations. The former are interpreted 
as high-order g-modes, the latter as low order p 
and/or g-modes. In the early nineties, the new 
OPA L opacity determinations ( Rogers fc Iglesiaj 
allowed to solve the long standing enigma 
of B type pulsators. Rosseland mean opacities 
computated from these tables were substancially 
increased in the Z-bump region compared to the 
previously computed opacities. This opacity en- 
hancement (as large as a factor of 3) provided a 
boost for the K-mechanism potentially activated 
in the metal opacity bump region and could then 
explain the excitation of pulsations in (3 Cephei 
and SPB theoretical standard models (ICox et al 



1992: Kiriakidis et all 1992: Moskalik et al 



19921: 



Dziembowski et al.l 119931^ ) However, as obser- 
vational capabilities progress, B type pulsators 
offer new challenges to the theorists. Some of 
them are no w found in low me t allicit y environ- 
ments (e.g. iKolaczkowski et "aL I l2006l and ref- 
erences therein). A large number of pulsation 
modes are indeed detected in B stars, with fre- 
quency spectra sometimes revealing new discrep- 
ancies between theory and observations. As an 
exa mple, the two Cephei stars, v Eri and 12 



J f^r — \^ 

Lac (jHandler et al.ll2006l : IJerzvkiewicz et al.ll2005l 



and references therein) present frequency spectra 
which cannot be reproduced wi th B type stars 



standard models (e.g. for i/ Eri, lAusseloos et al 



20041 ). iPamvatnvkh et all l|2004f) showed that an 



"ad hoc" factor 4 enhancement of the iron group 
elements, in the Z-bump region could help repro- 
ducing the observed pulsation spectrum of v Eri. 
The combined effects of radiative levitation and 
gravitional settling provides a natural explanation 
for this local metal enrichment and could then 
explain the observation of more than 60 /3 Cephei 
stars in low metallicity environments. 

Following the results of IPamvatnvkh et al 



(|2004l ). iBourge fc AlecianI (|2006l ) Dresented pre- 
liminary evolutionary models of /3 Cephei stars 
evolved from the ZAMS up to 1 Myr: their evolu- 
tionary models did not include consistent diffusion 
computations but they included, at each timestep, 
a rough estimation of the diffusion-induced iron 
accumulation in the Z-bump region. After one 
million years, their inferred iron enrichment is in- 
creased by a factor of 12 in the iron opacity bump. 
The stability analysis of these models showed that 
such iron enrichments makes the theoretical fre- 
quency spectra denser, which could help reducing 
the discrepancies between asteroseismic observa- 
tions and theoretical stability analysis. 

Following the same idea, iMiglio et all (|2007bl ) 
carried out a parametric study of the effects of lo- 
cal iron enhancement on the stability of SPB and 
(3 Cephei stars. The parametric iron accumulation 
profile introduced in their models is described by 
a gaussian function centered at logT ~ 5.2, cali- 
brated to be consi stent with the i r on acc umulation 
profile deduced bv lRichard et al. I (l200lh in their A 
and F star models. As Bourge fc AlecianI ( 20061 ). 
Miglio et all (|2007bf l showed that such accumu- 



lation profiles widen the theoretical SPB and /? 
Cephei instability strips and increases the number 
of excited modes. 

Bourge fc AlecianI (|2006h 



T he computat i ons by 



and iMigUo et al.1 (|2007br ) reduced the discrepan- 



cies between theoretical and observed pulsations. 
However they did not discuss the credibility of the 
iron profiles introduced in their models, nor the 
stability of the induced /i- gradients. 

On the other hand, Miglio et al. ( 2007a[ ) anal- 
ysed the effects of uncertainties in the opacity 
computations on the excitation of pulsation modes 
in B type stars. They coi npared models computed 
with OPAL opacities ( Iglesias fc Roger"^ 1996 ) 
with models comput ed with the re c ently updated 
OP opacity tables ( SeatonI 2005t Badnell et al 



20051 ). The OP updated data lead to an en- 



hancement o f 18% o f the opacity in the Z-bump. 



Miglio et al.l (|2007al ) showed that these differ- 



ences considerably affect the theoretical insta- 
bifity strips of /3 Cephei and SPB stars. Using 
the OP opacities drastically increases the num- 
ber of excited modes and significantly extends the 
theoretical instability strips. These new opacity 
computations could help solve, at least partly, 
the discrepancies between theoretical predictions 
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and observations in f3 Cephei, SPB and hybrid 
stars. Up to now, no evolutionary models have 
been computed for B-type stars including con- 
sistent diffusion computations and the new OP 
opaciti es. However in the light of iMiglio et al 



(l2007ah re sults, we expect th e iron accumulation 
derived bv lPamvatnvkh et al.l ( 2004h to recouncile 
observations and theory (by a factor of 4) to be 
overestimat ed. In this context t he iro n enrichment 
inferred by iBourge fc AlecianI (|2006l ) (by a factor 
of 12) for a 1 Myr-model should also be too large, 
which suggests the presence of a "forgotten" pro- 
cess (such as the thermohaline convection) which 
could reduce the diffusion-induced iron accumula- 
tion. 

5.2. The sdB stars 

Contrarily to SPB stars, subdwarf B stars (sdB 
stars) are evolved, compact objects with low-mass 
cores and outer H-rich envelopes. They are expe- 
riencing the core helium burning phase on the ex- 
tended horizontal branch. They all present chem- 
ical peculiarities attributed to gravitational set- 
tling and radiative Icvitation in the presence of 
weak stellar winds. Subdwarf B stars host two 
groups of pulsators: the first one, composed of 
hot sdB stars, is characterised by rapid ocillations, 
with periods ranging from 80 to 600s, caused by 
low-order, low-degree p-modes. The second group, 
composed of the coolest sdB stars, oscillates with 
periods in the 2000-9000s range, due to high-order, 
low-degree g-modcs. Both types of pulsators are 
driven by the K-mechanism acting in the iron-peak 
element opacity bump. 

The relation between the envelope metal con- 
tent and the puls ation driving m e chani sm was 
first established bv lCharpinet et aL for the 

short-period B subdwarf pulsators. They showed 
that a solar metal content in the envelope of their 
models is unable to excite modes but that mod- 
els with envelopes containing an overabundance of 
metals can globally destabilize modes such as ob- 
served. They argued that this metal enrichment 
needs only to occur in the driving region itself and 
not necessarily in the whole envelope as assumed 
in their models. 

A natural explanation for the required local 
metal enhancement is the radiative levitation pro- 
cess, which is expected to act efficiently in the ra- 
diative envelopes of subdwarf B stars. In these 



stars, as in A and F stars, the competition between 
gravitational settling and radiative aceleration is 
expected to produce local accumu lations of heavy 
elements. ICharpinet et al.l (jl997l ) presented more 
sophisticated and realistic models of sdB stars, 
the so-called "second generation models" , includ- 
ing non uniform iron abundances. These enve- 
lope models rely on the assumption that a state 
of diffusive equilibrium is reached between radia- 
tive levitation and gravitational settling operating 
on iron (disregarding other potentially competing 
proce sses) , leading to s t able n on uniform iron pro- 
files. ICharpinet et al. ( 1997t ) demonstrated that 
such local iron accumulations in the Z-bump re- 
gion may produce efficient excitation for pulsa- 
tion modes. Similar models were successfully 
used for detailed asteroseismic studies of several 
sdB stars: they globally nicely reproduced the 
global properties of the pulsating sdB stars and 
were able to closely reproduce their observed pe- 
riods (IBrassard et al.ll2001t ICharpinet et al.ll2001 



Charpinet erall (|2009l ). 



2005al ibl. l2006l . l2008l) 

As emphasized by 
these "second generation models" however suffer 
from several shortcomings. The theoretical insta- 
bility strip predicted by these models is larger than 
the observed one and the observed period ranges 
in individual stars are also usually narrower than 
predicted. These discrepancies are likely due to an 
overestimate of the iron abundance in the Z-bump 
region. 

In the "second generation models" of Charpinet 
et al., the opaci ties are computed using the OPAL 



opacity tables (jlglesias fc Roger 3 ll996l) . As dis- 
cussed previously for SPB and /? Ce phei stars, the 
use of the most rece nt OP opacities (jSeaton 120051: 
Badnell et al. 20051 ) may also have a direct im- 
pact on the properties of the sdB pulsation driv- 
i ng mechanism . Simi l arly to other types of stars 
(jjefferv fc Said 120071 : iMiglio et all l2007al ). using 
OP opacities instead of OPAL data, would in- 
crease the opacity in the Z-bump, resulting in a 
more efficient pulsation driving mechanism. 

Here again we find evidences that the iron ac- 
cumulation needed to account for the observations 
must be lower than obtained by the atomic diffu- 
sion alone. Introducing thermohaline convection 
could certainly help reducing the discrepancies. 
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5.3. The 7 Doradus stars 

The 7 Doradus stars are A-F main sequence 
pulsators. Their osciUations arc interpreted as 
high-order g-modes with periods ranging from 0.35 
to 3 days. They are driven by a flux blocking 
mechanism at the base of the i r conv ective enve- 
lope, as shown bv lGuzik et al. ( 20001) with fr ozen 
convection models a nd bv lDupret et al.l (j2004t ) and 
Dupret et al. I (l2005h with time-dependent convec- 
tion models. It has been demonstrated that the 
"convection blocking" of radiation can drive high 
order g-modes only in stellar models for which the 
temperature at the bottom of the convective enve- 
lope, ranges between 2 x lO^K (which corresponds 
nearly to the Z-bump region central temperature) 
and 4.8 x lO^K. The depth of the convective enve- 
lope plays a crucial role in the d riving mechanism 
proposed by Guzik et al. ( 200G| ) to explain the 7 
Doradus pulsations. 

As a consequence, the seismic properties of 7 
Doradus stars are quite sensitive to atomic diffu- 
sion including radiative levitation because of the 
way they in fluence the convective z ones. As em- 
phasized bv lMontalban et al.l (j2007l ). the diffusion 
effects on convection are of various kinds. The He 
settling which leads to H-enrichment in the exter- 
nal stellar layers induces an opacity increase which 
affects the depth of the surface convective zone. 
On the other hand, as discussed in previous para- 
graphs, radiative levitation produces reservoirs of 
iron-peak elements which may induce the appear- 
ance of an iron convective zone near 2 x lO^K. 
Consequently the driving mechanism of 7 Doradus 
stars is expected to be sensitive to the thermoha- 
line convection which occurs below the heavy ele- 
ments accumulation zone. 

6. Conclusion 

In the present paper, we tested the effect of 
thermohaline convection on the amplitude of iron 
accumulation which occurs inside A-F type stars 
due to atomic diffusion. We computed two se- 
ries of models for three stellar masses: 1.5, 1.7 
and 1.9M0. In both series, atomic diffusion was 
consistently computed for the elements H, He, C, 
N, O, Ca and Fe. The computations included 
the radiative accelerations on C, N, O, Ca and 
Fe, by us ing the semi-analyt i cal pr escription pro- 
posed bv lAlecian fc LeBland (|2002l ) and improved 



bv lLeBlanc fc AlecianI (|2004l ). 

No thermohaline convection was introduced 
in the first series of models. We verified that 
our results w ere s imilar to those obta i ned b y 
Richer et al.l (|2000l ) and iRichard et all (|200ll ). 



with a large accumulation of iron in the iron peak 
elements opacity bump region. Thermohaline con- 
vection was then introduced for the second series 
of models. Considering the short time scales of 
this mixing process compared to the evolution 
time scales, and the fact that mixing occurs until 
the unstable ^-gradients disappear, we introduced 
it in such a way that the //-gradients were kept 
close to zero or negative all along the evolutionary 
tracks. In practice, we used an extra diffusion co - 
efficicnt, as proposed by Kippenhahn et al. ( 198G| ). 
which is large enough to reach this purpose. 

We presented in detail the results obtained for 
the 1.7Mq case, during the main sequence phase. 
We also showed some results for the two other 
stellar masses considered here. The effect of ther- 
mohaline convection is very efficient in strongly 
reducing the heavy element accumulation in the 
opacity bump region. Whereas the iron abundance 
can increase by up to factors larger than 100 when 
thermohaline convection is neglected, it never in- 
creases by more than a factor of 15 when it is taken 
into account. 

It is important to note that some iron accumu- 
lation always remains, even in the presence of ther- 
mohaline convection. The reason why this hap- 
pens is directly related to helium settling, which 
leads to a stabilizing /(-gradient. This stabilizing 
effect allows in turn some heavy material abun- 
dance increase up to the point where the /i-profile 
becomes flat. 

As a consequence, the cases where iron accumu- 
lation has been invoked to account for speciflc ob- 
servations in stars remain. Indeed, in many cases 
the observations would be better explained with 
a smaller iron accumulation than that given by 
atomic diffusion alone. All these cases (Am stars, 
SPB stas, sdB stars, /? Cephei stars and 7 Do- 
radus stars as discussed in Section 5) will have to 
be studied individually for a detailed comparison 
with the observations. 

In the present computations, nickel was not in- 
cluded. Although less important than iron, it may 
signiflcantly contribute to the opacity and will 
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have to be included in future work. No mass loss 
nor extra turbulence were taken into account. In- 
cluding such macroscopic motion may reduce the 
heavy element accumulation even without thermo- 
haline convection. However, thermohaline convec- 
tion is a physical process which must occur in any 
case every time the /^-gradient is unstable. The 
final accumulation of elements thus strongly de- 
pends on this process, even in the presence of other 
macroscopic motion. 

In this paper, we have demonstrated the im- 
portance of thermohaline convection which has to 
be taken into account in the computation of the 
abundance variations induced by atomic diffusion. 
Precise applications to specific stellar cases still 
have to be undertaken. They will lead to im- 
portant consequences on the abundance anomalies 
observed in stars, as well as their oscillating prop- 
erties through several processes like K-mechanism, 
convective blocking, etc. The possible applications 
of thermohaline convection in stellar astrophysics 
are vast and will certainly be a topic of interest in 
years to come. 

REFERENCES 

Alecian, G., & LeBlanc, F. 2002, MNRAS, 332, 
891 

Ausseloos, M., Scufiaire, R., Thoul, A., & Aerts, 
C. 2004, MNRAS, 355, 352 

BadncU, N. R., Bautista, M. A., Butler, K., De- 
lahaye, F., Mendoza, C, Palmeri, P., Zeippen, 

C. J., & Seaton, M. J. 2005, MNRAS, 360, 458 

Balmforth, N. J., Cunha, M. S., Dolcz, N., Gough, 

D. O., & Vauclair, S. 2001, MNRAS, 323, 362 

Bourge, P.-O., & Alecian, G. 2006, ASP Confer- 
ence Series, vol. 349, 201 

Brassard, P., Billcres, M., Charpinet S., Liebert, 
J., & Saffer, R. A. 2001, ApJ, 563, L1013 

Charbonnel, C., Zahn, J. P., 2007, A&A, 476, L29 

Charpinet S., Fontaine, G., Brassard, P., & Dor- 
man, B. 1996, ApJ, 471, L103 

Charpinet S., Fontaine, G., Brassard, P., Chayer, 
P., Rogers, F. J., Iglesias, C. A., & Dorman, B. 
1997, ApJ, 483, L123 



Charpinet, S., Fontaine, G., & Brassard, P. 2001, 
PASP, 113, 775 

Charpinet S., Fontaine, G., Brassard, P., Green, 
E. M., & Chayer, P. 2005a, A&A, 437, 575 

Charpinet S., Fontaine, G., Brassard, P., Billcres, 
M., Green, E. M., & Chayer, P. 2005b, A&A, 
443, 251 

Charpinet S., Silvotti, R., Bonanno, A., Fontaine, 
G., Brassard, P., Chayer, P., Green, E. M., 
Bergeron, P., et al. 2006, A&A, 459, 565 

Charpinet S., Van Grootel, V., Reese, D., 
Fontaine, G., Green, E. M., Brassard, P., & 
Chayer, P. 2008, A&A, 489, 377 

Charpinet S., Fontaine, G., & Brassard, P. 2009, 
A&A, 493, 595 

Cox, A. N., Morgan, Siobahn, M., Rogers, F. J., 
Iglesias, & C. A. 1992, ApJ, 393 272 

Dupret, M.-A., Grigahccnc, A., Garrido, R., 
Gabriel, M., & Scufiaire, R. 2004, A&A, 414, 
L17 

Dupret, M.-A., Grigahcene, A., Garrido, R., 
Gabriel, M., & Scufiaire, R. 2004, A&A, 435, 
927 

Dziembowski, W. A., & Pamyatnykh, A. A. 1993a, 
MNRAS, 262, 204 

Dziembowski, W. A., Moskalik, P., & Pamyat- 
nykh, A. A. 1993b, MNRAS, 265 588 

Gargett, A., Ruddick, B., eds. 2003, Double- 
Diffusion in Oceanography (Oxford:Pergamon), 
381 

Grevesse, N., & Noels, A. 1993, in Origin and Evo- 
lution of the Elements, Eds Kubono, S. and Ka- 
jino, T., 14 

Guzik, J. A., Kaye, A. B., Bradley, P. A., Cox, A. 
N., & Neuforge, C. 2000, ApJ, 542, L57 

Handler, G., Jerzykiewicz, M., Rodrguez, E., Uyt- 
terhoeven, K., Amado, P. J., Dorokhova, T. 
N., Dorokhov, N. I., Poretti, E., Sareyan, J. -P., 
Parrao, L., and 21 coauthors, 2006, MNRAS, 
365, 327 

Hui-Bon-Hoa, A. 2008, Ap&SS, 316, 55 



18 



Iglesias, C. A., & Rogers, F. J. 1996, ApJ, 464, 
943 

Jeffery, C. S., & Saio, H. 2007, MNRAS, 378, 379 

Jerzykiewicz, M., Handler, G., Shobbrook, R. R., 
Pigulski, A., Medupe, R., Mokgwetsi, T., Tl- 
hagwane. P., & Rodrguez, E. 2005, MNRAS, 
360, 619 

Kippenhahn, R., Ruschenplatt, G., & Thomas, 
H.-C. 1980, A&A, 91, 175 

Kiriakidis, M., El Eid, M. F., & Glatzel, W. 1992, 
MNRAS, 255, 1 

Kolaczkowski, Z., Pigulski, A., Soszyiiski, I., Udal- 
ski. A., Kubiak, M., Szymanski, M., Zebrun, 
K., Pietrzyhski, G., Wozniak, P. R., Szewczyk, 
O., & Wyrzykowski, L. 2006, Memorie della So- 
cieta Astronomica Italiana, vol. 77, 336 

LeBlanc, F., & Alccian, G. 2004, MNRAS, 352, 
1329 

Michaud, G. 1970, ApJ, 160, 641 

Miglio, A., Montalban, J., & Dupret, M.-A. 2007, 
MNRAS, 375, 21 

Miglio, A., Bourge, P.-O., Montalban, J., & 
Dupret, M.-A. 2007, CoAst, 150, 209 

Montalban, J., Miglio, A., & Theado, S. 2007, 
CoAst, 150, 137 

Moskalik, P., & Dziembowski, W. A. 1992, A&A, 
256, L5 

O'Donoghue, D., Lynas-Gray, A. E., Kilkenny, D., 
Stobie, R. S., & Koen, C. 1997, MNRAS, 285, 
657 

Pamyatnykh, A. A., Handler, G., & Dziembowski, 
W. A. 2004, MNRAS, 350, 1022 

Richard, O., Michaud, G., & Richer, J. 2001, ApJ, 
558, 377 

Richard, O., Theado, S., & Vauclair, S. 2004, 
SoPh, 220, 243 

Richer, J., Michaud, G., & Turcotte, S. 2000, ApJ, 
529, 338 

Rogers, F. J., & Iglesias, C. A. 1992, ApJS, 79, 
507 



Santos, N.C., Israelian, G., & Mayor, M. 2001, 
A&A, 373, 1019 

Seaton, M. J. 2005, MNRAS, 362, LI 

Stancliffe, R.J. 2009, MNRAS, 394, 1051 

Stancliffe, R.J., Glebbeek, E., Izzard, R.G., & 
Pols, O.R. 2007, A&A, 464, 57 

Stanchffe, R.J., & Glebbeek, E. 2008, MNRAS, 
389, 1828 

Stankov, A., & Handler, G. 2005, ApJS, 158, 193 

Stothers, R., & Simon, N.R. 1969, ApJ, 157, 673 

Theado, S, Alecian, G., & LeBlanc, F. 2009, in 
prep. 

Thompson, LB., Ivans, I.I, Bisterzo, S., Sneden, 
C., Gallino, R., Vauclair, S., Burley, G.S., 
Shectman, S.A., & Preston, G.W. 2008, ApJ, 
677, 556 

Turcotte, S., Richer, J., Michaud, G., & 
Christensen-Dalsgaard, J. 2000, A&A, 360, 603 

Ulrich, R. K. 1972, ApJ, 172, 165 

Vauclair, S. 1975, A&A, 45, 233 

Vauclair, S. 2004, ApJ, 605, 874 



This 2-column preprint was prepared with the AAS IKT^jpi 
macros v5.2. 



19 



